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Abstract. The first-year WMAP data taken at their face value hint that the Universe 
might be slightly positively curved and therefore necessarily finite, since all spherical 
(Clifford-Klein) space forms M 3 = S 3 /T, given by the quotient of <S 3 by a group 
r of covering transformations, possess this property. We examine the anisotropy of 
the cosmic microwave background (CMB) for all typical groups T corresponding to 
homogeneous universes. The CMB angular power spectrum and the temperature 
correlation function are computed for the homogeneous spaces as a function of the 
total energy density parameter Sltot in the large range [1.01, 1.20] and are compared 
with the WMAP data. We find that out of the infinitely many homogeneous spaces 
only the three corresponding to the binary dihedral group T* , the binary octahedral 
group O*, and the binary icosahedral group I* are in agreement with the WMAP 
observations. Furthermore, if fitot is restricted to the interval [1.00,1.04], the space 
described by T* is excluded since it requires a value of Sltot which is probably too 
large being in the range [1.06, 1.07]. We thus conclude that there remain only the two 
homogeneous spherical spaces S 3 fO* and S 3 / 1* with f^tot of about 1.038 and 1.018, 
respectively, as possible topologies for our Universe. 



PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 



1. Introduction 



At present, all data are consistent with, and in fact strongly support, the standard big- 
bang model in which the time evolution of the Universe is described by the Friedmann- 
Lemaitre-Robertson- Walker metric. Accordingly, the Universe possesses the space-time 
structure R x M where M. describes the "space" of cosmic time, and Ai the three- 
dimensional comoving space section of constant curvature K = +1,0 and —1. The 
Einstein field equations respectively the Friedmann equations for the cosmic scale factor 
do not fix the curvature a priori. Instead, the curvature parameter K has to be inferred 
from a determination of the total energy density e to t of the Universe via the relation 
(c = 1) K = H$al(Cl tot - 1), fi to t := £tot/ £ crit, where e CTit := |^ denotes the critical 
energy density, a the cosmic scale factor, and H the Hubble constant (all quantities 
at the present epoch). Furthermore, it is a mathematical fact, although not always 
appreciated (see, however, the remark on early works below), that fixing the curvature 
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K does not determine uniquely the global geometry of A4, i.e. the topology and thus 
the shape of the Universe. Only if it is assumed that the Universe is simply-connected, 
the possible homogeneous 3-spaces Ai of constant curvature K are given by the 3-sphere 
S 3 (K = +1), Euclidean 3-space S 3 (K = 0), or hyperbolic 3-space H 3 (K = — 1). In this 
case, the Universe is finite for positive curvature (fitot > 1) and infinite for vanishing 
(fi tot = 1) or negative curvature (fi tot < 1). However, most 3-spaces M. of constant 
curvature are multi-connected and are given by the quotient of S 3 , S 3 , or 7i 3 by a group 
T of covering transformations, i.e. M. = S 3 /T, S 3 /T, or 7i 3 /V. In this case, the Universe 
is again finite for positive curvature, but can be finite too if it is flat or negatively curved. 

Here, we would like to remark that the question whether the space of the Universe 
is finite and possibly multi-connected has been discussed during the last century by 
several cosmologists, e.g. by Schwarzschild pQ, Einstein j2j, Friedmann 011]) Lemaitre 
Heckmann and Schucking and Ellis [Zj, to mention only a few. 

The concordance model of cosmology (ACDM model) assumes a flat Universe 
with the topology of £ 3 with a positive cosmological constant A, i.e. £Ia '■= t^W = 
1 — f2 ma t ~~ ^rad with f2 mat = f2bar + ^cdm, where the various ^-parameters denote the 
present value of the baryonic (bar), cold dark matter (cdm), matter (mat) and radiation 
(rad) energy densities in units of e cr i t . Three variants of the concordance model have 
been presented by the WMAP team jHj providing a good overall fit to the temperature 
fluctuations 5T of the cosmic microwave background radiation (CMB) on small and 
medium scales, but there remains a strange discrepancy at large scales as first observed 
by COBE [H| and later substantiated by WMAP PU|. 

The suppression of the CMB anisotropy at large scales respectively low multipoles 
can be explained if the Universe is finite. Recent analyses concerning the suppression 
at low multipoles in the WMAP data can be found in [TT1 H2J EH EU EH] . As discussed 
before, a finite Universe is naturally obtained, if the total energy density exceeds the 
critical value one, i.e. Q t ot > 1- Interestingly enough, the WMAP team reported 
PH] fitot = 1.02 ± 0.02 together with Q har = 0.044 ± 0.004, fi mat = 0.27 ± 0.04, and 
h = 0.7llo!o3 f° r ^ e P resen t day reduced Hubble constant (the errors give the lcr- 
deviation uncertainties only). Taking at their face value, these parameters hint to a 
positively curved Universe possessing the geometry of S 3 or of one of the spatial space 
forms iS 3 /r. (One should keep in mind, however, that the WMAP values depend on 
certain priors, and, furthermore, include the la-errors only. Thus it would be too early 
to conclude that the data definitively exclude a negatively curved Universe. In fact, 
we have recently shown ^21 EH] that the non-compact, but finite hyperbolic Picard 
universe describes well the CMB anisotropy and the observed suppression of power at 
large scales.) 

In this paper, we present a systematic comparison of the predictions with the CMB 
anisotropy for universes possessing homogeneous spherical topology. This comparison is 
made possible for two reasons. First of all, the spherical spaces were classified already by 
1932 fOlEI and thus their mathematical structure is known. (This is in contrast to the 
case of hyperbolic manifolds which are not yet completely classified; even the manifold 
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with the smallest volume is not yet known.) Second, due to an efficient numerical 
algorithm described in our recent paper ^3], we are able to take in the Sachs- Wolfe 
formula a large number of vibrational modes into account and thus to predict sufficiently 
many CMB multipoles for the various spherical spaces which in turn allow a detailed 
comparison with the WMAP data. Since the CMB spectrum depends sensitively on 
the curvature radius, the comparison is performed as a function of Q tot hi the large 
interval [1.01, 1.20] in order to determine for a given spherical space form the best- 
fitting value of the total energy density - under the condition, of course, that the space 
under consideration is able to describe the data at all. 

Recently, Luminet et al. ^T] proposed the Poincare dodecahedron, which is one 
of the well-known spherical space forms (see section |21 for details), as a model for the 
geometry of the Universe. In their preliminary study involving only the three lowest 
multipoles (/ = 2,3,4) they found, indeed, for f2 tot = 1.013 a strong suppression of the 
CMB power at I = 2 and a weak suppression at I = 3 in agreement with the WMAP 
data. However, in jTI] only the first three vibrational modes of the dodecahedral space 
with wave number f3 = 13, 21 and 25 (comprising in total 59 eigenfunctions) have been 
used, and there thus remained the question about how this extremely low wave number 
cut-off affects the predictions of the multipoles, since experience shows that increasing 
the cut-off usually enhances the integrated Sachs- Wolfe and Doppler contributions. 

In our recent paper ^3] we presented a thorough discussion of the CMB anisotropy 
for the dodecahedral space topology using the first 10521 eigenfunctions corresponding to 
the large wave number cut-off f3 = 155. The contributions of higher wave numbers up to 
(3 = 1501 were taken into account with respect to their mean behaviour. Taking within 
the tight-coupling approximation not only the ordinary, but also the integrated Sachs- 
Wolfe and also the Doppler contribution into account, we were able to predict sufficiently 
many multipole moments such that a detailed comparison of the dodecahedral space 
model with the WMAP data could be performed. We found that the temperature 
correlation function for the dodecahedral universe possesses very weak correlations at 
large scales in nice agreement with the WMAP data for Q tot in the range 1.016 . . . 1.020. 
There thus arises the interesting question whether the dodecahedral space is the only 
spherical space form able to describe the CMB data. In ^H] the CMB anisotropy is also 
studied for the dodecahedron, the binary octahedral group, and the binary tetrahedral 
group. 

The main result of the present paper is that while almost all homogeneous spherical 
spaces have to be excluded as possible geometries for the Universe, there is one particular 
space form, defined by the binary octahedral group, which agrees for Q t ot — 1-038 
with the WMAP data even better than the dodecahedron. We observe that the best- 
fitting values obtained for f2 tot are different for the binary octahedral space and the 
dodecahedron, but both values lie well within the la-band determined by WMAP. It 
remains to be seen whether future data will enable us to definitely eliminate one of the 
two space forms in favour of the other as describing the true topology of the Universe. 

Our paper is organised as follows. In section |2l we summarize the main properties 
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of the existing homogeneous spherical space forms and of their vibrational modes. Our 
main results are presented in section [3] which contains a detailed comparison with the 
WMAP data of the CMB angular power spectrum for the various types of spherical 
spaces. In addition to the power spectrum, we study also the so-called S(p) statistic 
[TU] which measures the suppression at large angular scales directly in terms of the 
temperature correlation function. Section 0] contains our conclusions. 

2. The spherical space forms and their vibrational modes 

In section 2 of ^1] we have already described the three-dimensional spaces M. of constant 
positive curvature K = 1, and therefore we refer the reader to this paper for details. 
The spherical spaces were classified by 1932 [TBI E! an d are given by the quotient 
M. = S 3 /T of the three-sphere S 3 under the action of a discrete fixed-point free subgroup 
T C SO (4) of the isometries of S 3 . All these manifolds are compact possessing the 
volume V(S 3 /T) = V(S 3 )/N, where N is the order of the group T, and are, apart from 
the universal covering space S 3 with volume V(S 3 ) = 27r 2 , multi-connected. To define 
the discrete fixed-point free subgroups T C SO (4) of isometries of S 3 , one makes use 
of the fact that the unit 3-sphere S 3 can be identified with the multiplicative group of 
unit quaternions {q}. The latter are defined by q := w + xi + y] + zi}, (w, x, y, z) G M 4 , 
having unit norm, \q\ 2 = w 2 + x 2 + y 2 + z 2 = 1. Here, the 4 basic quaternions {1, i, j, ij} 
satisfy the multiplication rules i 2 = j 2 = — 1 and ij = — ji plus the property that i and j 
commute with every real number. The distance d(qi,q 2 ) between two points q\ and q2 
on S 3 is given by cosc^, q 2 ) = WiW 2 + x x x 2 + y\y2 + Z\z 2 . 

The group SO(4) is isomorphic to S 3 x <S 3 /{±(1, 1)}, the two factors corresponding 
to the left and right group actions. In this paper, we are only interested in homogeneous 
manifolds M. = S 3 /T, in which case the group T contains only right-handed Clifford 
translations 7 G T that act on an arbitrary unit quaternion q G S 3 by left-multiplication, 
q — > jq, and translate all points q±,q 2 G S 3 by the same distance Xi l - e - ^(<?i;7<?i) = 
d(q 2 , , yq 2 ) = \. The right-handed Clifford translations act as right-handed cork screw 
fixed-point free rotations of S 3 . The following groups lead to homogeneous manifolds 
M = s 3 /r [na [12 [TBI [EH]: 

• The cyclic groups Z m of order m (m > 1). 

• The binary dihedral groups D\ m of order 4m (m > 2). 

• The binary tetrahedral group T* of order 24. 

• The binary octahedral group O* of order 48. 

• The binary icosahedral group I* of order 120. 

In tabled we give the right-handed Clifford translations which generate the above groups 

r. 

The vibrations on the homogeneous spherical spaces M. = S 3 /T are determined by 
the regular solutions of the Helmholtz equation 

(A + 0^(g) = , qeM , (1) 
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Table 1. The generators 71 and 72 for the groups r (a = (1 + V5)/2). 



satisfying the fundamental periodicity conditions 

^'\lkq) = <' J (g) , Vg G A* , V 7fc G T . (2) 

Here A denotes the Laplace-Beltrami operator on S 3 . The eigenfunctions on M. satisfy 
the orthonormality relation 

/ d^if)^\q)if)^{q) = -j- / d l i^{q)^{q) = Spp, Sii> . (3) 

The spectrum on M. is discrete, and the eigenvalues can be expressed in terms of the 
wave number (3 G N as Ep = (3 2 — 1 and are independent of the degeneracy index 
i = 1, . . . ,r M (/3), where r M (fi) denotes the multiplicity of the mode (3. It should be 
noted that for a given manifold M. the wave numbers /3 do not take all values in N. The 
allowed (3 values together with their multiplicities r M (f3) are explicitly known 1201121], 
see Table 121 

The eigenfunctions if) a ,l (q) on M. can be expanded into the eigenfunctions 4>pi m (q) 
on S 3 

v#%) = EEWtt • (4) 

;=0 m=-l 

Since the eigenfunctions on S 3 are explicitly known [2211211121112312111235 it remains to 
determine the expansion coefficients £,pi m (M) which satisfy as a consequence of eq. © 
the normalization condition 

E E (W^))* = AT . (5) 

1=0 m=-l 

In [T3] we have described our numerical algorithm to compute the coefficients £L m (.M). 
It uses a collocation method by imposing the periodicity condition (J2J. Using this 
method, we have computed in ^1] the expansion coefficients for T = I*, i. e. for the 
Poincare dodecahedral space T> = S 3 /I*, for (3 < 155 comprising the first 10521 
eigenfunctions. In addition, we have computed in [Hj the coefficients £gj m (.M) for 
/5 < 33 for some cyclic groups, some binary dihedral groups, the binary tetrahedral 
group, and the binary octahedral group. Based on these numerical results, we stated 
for homogeneous space forms the following relation as a conjecture (0 < I < f3 — 1) 

^ y, (^j.m)T = n t w , (6) 

^ m=-l i=l KHJ 
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where r s (f3) = /3 2 denotes the multiplicity of the vibrational modes on S 3 . The relation 
(0) has been found to hold within a numerical accuracy of 13 digits. However, for the 
inhomogeneous lens spaces [2H] L(12,5) and L(72,17) we have found that the relation 
(jHJ) does not hold. We thus concluded in [TJj that the relation (JHJ) is only valid for 
homogeneous 3-spaces. Recently, Gundermann ^H] provided a proof of our conjecture 
©• 

In the next section, we shall use the eigenfunctions (@J) to calculate CMB sky maps 
and the variance of the CMB anisotropy and relation (JHJ) to calculate the mean value of 
the CMB anisotropy for a variety of spherical spaces. 



3. The angular power spectrum ST 2 and the correlation function C(d) for 
spherical spaces 

The relative temperature fluctuations ST ^ of the CMB are caused by several effects 
which we shall compute within the tight-coupling approximation along the lines 
described in detail in Section 2 of [Ej. The dominant contribution at large scales 
is given by the ordinary Sachs-Wolfe (SW) effect which is a combination of the 
gravitational potential <&(rj,T,6,(j)) at the surface of last scattering (SLS), and the 
intrinsic temperature fluctuation |<5 7 (r?, r, 9, 0) due to the imposed entropic initial 
conditions, where <5 7 denotes the relative perturbation in the radiation component. (Here 
r) denotes the conformal time and (r(rj),0,<f)) the spherical coordinates of the photon 
path in the direction n = n(8, <ft), where we assume that the observer is at the origin of 
the coordinate system, i.e. at (r b s , 6* bs, 0obs) = (0,0,0).) 

The gravitational potential $ is identified with the (scalar) perturbation of the 
Friedmann-Lemaitre-Robertson- Walker metric which for an energy-momentum tensor 
T^v with Tij = for i ^ j (i,j = 1, 2, 3) can in conformal Newtonian gauge be written 
as [HI 

ds 2 = a 2 (7]) [(l + 2$)dr/ 2 -(l-2$)|rff| 2 ] . (7) 

Here a(t]) denotes the cosmic scale factor as a function of conformal time t] and \dx\ 2 
the line element on S 3 

\dx\ 2 = dr 2 + sin 2 r (d6 2 + sin 2 9 dcf) 2 ) (8) 

with 0<r<7r, 0<6 <n, O<0<2vr. 

The ordinary Sachs- Wolfe (SW) contribution to the temperature fluctuation is given 

by 

5T sw (n) 1 

^ = $(??SLS,^SLS,#,0) + ^(^SLS^SLS,^) (9) 

with tsls : = Vo~ Vsls, where rj and tjsls denote the conformal time at the present epoch 
and at the time of recombination corresponding to a redshift ^sls = 1089, respectively. 
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For a given spherical space Ai, the metric perturbation can be written as an expansion 
in the eigenf unctions on M. 

r M (P) 

*fo,r,M) = E' E , (io) 

/3>3 i=l 

where in the mode summation only the modes with (3 > 3 have been taken into account 
(if they exist), since the wave numbers (3 = 1,2 correspond to modes which are pure 
gauge terms [2E]. The prime in the summation over the modes f3 indicates that the 
spectrum of a given manifold Ai does not contain all /3 G N, see Table El The functions 
$p(T]) determine the time evolution and will be factorized $^(r/) = 3>]g(0) 573(77) with 
5/3(0) = 1. The functions 573(77) do not depend on the degeneracy index i, since the 
associated differential equation depends only on the eigenvalue E^ 4 which is independent 
of i. The initial values $«(0) are the primordial fluctuation amplitudes and are assumed 
to be Gaussian random variables with zero expectation value and covariance 



^(0)^,(0)^ = 8pp8«P*{p) . (11) 

Here -P$(/3) denotes the primordial power spectrum that determines the weight by which 
the primordial modes (3 are excited, on average. The average (...) in (fTTj) denotes an 
ensemble average over the primordial perturbations which are supposed to arise from 
quantum fluctuations, by which the Universe is "created". In the following, we shall 
assume that the primordial power spectrum is in good approximation described by the 
scale-invariant Harrison-Zel'dovich spectrum 

p M - J^T) ■ < 12 > 

Here a is a normalization factor which will be determined from the CMB data. 

The temperature fluctuations ST(n) of the microwave sky can be expanded into real 
spherical harmonics Yi m (n) on S 2 , 

00 1 

ST(h) := E E a i™Yirn(h) , (13) 

1=2 m=-l 

where the monopole and dipole terms, I = 0, 1, are not included in the sum (|13j) . From 
the real expansion coefficients a\ m one forms the multipole moments 

c '^(tw') a*) 

m=— I 

and the angular power spectrum 

8Tf := ^±12 Cl . (15) 

The average (...) in (|14j) denotes an ensemble average over the primordial perturbations 
as in eq. (jll|) . respectively an ensemble average over the universal observers. 
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one arrives with the help of relation (jUJ) at the following expression for the 
ordinary Sachs- Wolfe contribution to the multipole moments for a given spherical space 



Here Rpi(r) denote the "radial functions" on S 3 which can be expressed in terms of 
Gegenbauer polynomials, see eq.(10) of [Hj. The expression (|T7|) shows in a transparent 
way, why the lowest multipoles are in general suppressed for the multi-connected 
spherical spaces, the more the more wave numbers (3 are missing in the vibrational 
spectrum. Let us consider, for example, the quadrupole moment, I = 2. Then the 
summation over modes in (|17j) runs for the simply-connected manifold S 3 over all natural 
numbers with (3 > 3. In contrast, for the Poincare dodecahedral space T> there is a large 
gap between (3 = 3 and 12, since the lowest contributing mode occurs only at (3 = 13, and 
thus the missing modes lead to a suppression. The suppression on T> gets even stronger, 
because there exist only the three modes (3 = 21, 25 and 31 in the wave number interval 
13 < (3 < 33 (see table EJ). 

Now, we would like to discuss this suppression mechanism in more detail. Tabled 
shows that the allowed wave numbers (3 are for all homogeneous spherical spaces given 
by odd natural numbers, except for the homogeneous lens spaces S 3 /Z m with m odd 
> 1 for which (3 runs through all natural numbers with (3 > m + 1, in addition to the 
lowest odd wave numbers between 1 and m. While the (3 spectrum consists for the 
homogeneous lens spaces Z m with m even > 2 of all odd natural numbers, there are for 
all other spherical spaces at low wave numbers below a given threshold (3 t h "gaps" in the 
spectrum. As can be seen from Table El the threshold value (3^ increases if the volume 
V{M) decreases, i. e. (3 t h = 4[(m + l)/2] + 1, 13, 25 and 61 for the spaces belonging to 
the groups D\ m , T*, O* and I*, respectively. This demonstrates clearly the important 
role played by the spatial volume, i. e. the topology of the Universe. 

In addition to the wave number gaps, there is another important imprint of topology 
on the multipole moments C\ due to the multiplicities r M (f3), also given in Tabled by 
which the vibrational modes are weighted in the mode sum (|17j) . Again, this effect is 
significant for the low vibrational modes and thus for the low multipole moments which 
can be considered as carrying the fingerprints of the topology of the Universe. 

On the other hand, for the large multipole moments, I 3> (3 t h, the details of the 
different topologies get washed out due to the identical mean asymptotic behaviour of 
the multiplicities for all multi-connected spaces M. (except Z m , m odd) given by 



M = S 3 /T 




(17) 



V{M) 



(3 2 + . . . for (3 



(18) 



oo 
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Table 2. The eigenvalue spectrum for the groups T. 



which leads to the universal formula (7 3> (3th) 

2 
9 



9 

Cf^) * q E ^(/3)^(%Ls)^(r S Ls) • (19) 

/3odd,/3>/ 



For the homogeneous lens spaces Z m , m odd > 1, one obtains instead 

r M {(3) = 1/5 2 + ... for - oo , (20) 
and thus from (JT7j) for Z 3> /3th 

1 oo 

Q sw (Z m ) ~ - MP) qKvsls) RU^s) (21) 

/3=«+l 

which should be numerically not very different from the expression ()19|1 where the 
additional factor 2 should recompensate for the missing even /3-values in the sum f|19|) . 
The asymptotic behaviour given in eqs. (jl8j) and (f2T)J) is in agreement with Weyl's law, 
as can be seen as follows (k — > oo): 



M M {k) :=#{(3\P<k} = J^r M ([3) 



/3<k 



v iV / 6?r 

/3odd,/3<fc v 7 

for T Z m , m odd > 1, and 

^ m w = e0 2 ) + - = iP* 3 + '- (23) 

(3=1 V ' 

for T = Z m ,m odd > 1. One thus sees that the ordinary Sachs- Wolfe contribution to 
the large multipoles is identical for all spherical spaces, including S 3 , in agreement with 
the expectation that the cosmic topology is most clearly seen at large scales. 

The above discussion showed how the topology of the Universe influences via the 
vibrational modes the CMB anisotropy at large scales. The whole story is, however, 
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more subtle, since an additional /-dependence comes in eq.(JT7J) from the radial functions 
R/3i(rshs)i an d furthermore, there is an important dependence of the multipoles on the 
curvature radius respectively on Q tot which is determined by both the time evolution via 
^(^sls) an d the radial functions. The interplay of all these effects is responsible for the 
rather complicated structure displayed by the numerical computations to be illustrated 
below. 




1.15 n tot l 



Figure 1. The binary tetrahedral group T*. Panel (a) shows the f2tot dependence 
of the mean value of the first three angular power moments <5T ; 2 . The horizontal 
lines indicate the corresponding WMAP values for 5Tf. Panel (b) displays the S(p) 
statistics (p = 60° full curve and p = 20° dotted curve) in units of pK 4 in dependence 
on Otot' The corresponding WMAP values are indicated as horizontal lines. 
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Figure 2. The same as in figure ^ for the binary octahedral group O* . 



Let us now discuss the results for the different spherical space forms. In all these 
computations the CMB anisotropy is obtained using the complete Sachs- Wolfe formula 
{r{ri) := r] - 77) 



T y 



)=£' E 

/3>3 i=l 



dr 
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Figure 4. The same as in figureHfor the binary dihedral group D%, i.e. m = 2. 
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Figure 5. The same as in figure \l\ior the binary dihedral group D^q, i. e. m = 5 



f (/9) 



+ 2 E £ / *iT**' H * MI • (24) 

/3> 3 i=l J, ?SLS ' 

The details of the computation of the quantities needed in are described in |12J. 
The first two terms in (J24j) are the ordinary Sachs- Wolfe contribution Q discussed 
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Figure 6. The same as in figure ^ for the cyclic group Z\, i. e. for the space S 3 . 
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Figure 7. The same as in figure ^ for the cyclic group Z2, i.e. for the projective 
space V 3 ■ 
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Figure 8. The same as in figure ^ for the cyclic group Zq. 
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above. The next term involving the spatial covariant divergence of the velocity field 
is the Doppler contribution. The integral over the photon path yields the integrated 
Sachs- Wolfe contribution. The /3 max cut-off is chosen sufficiently high in order to obtain 
enough multipoles which then enables us to normalise the 5Tf spectrum in the range 
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/ = 20 to 45 according to the WMAP first-year data [TO]. In our computations we use 
Anax = 3001 for O tot < 1.05, /3 max = 2001 for 1.05 < Q tot < 1.1, and /3 max = 1501 for 
1.1 < Q t ot < 1.2. In the following we fix the cosmological parameters as f2b ar = 0.046, 
^mat — 0.28, and h = 70, in agreement with the WMAP data. The free parameter is 
Qa = ^tot — ^mat — ^rad, i- e. the density corresponding to the cosmological constant. 
Although f^A is varied, we show in the following figures f2 tot on the abscissa. 

In figures [U to [5] we show the dependence of the large scale CMB anisotropy on 
fl to t for the binary tetrahedral space, the binary octahedral space, the dodecahedral 
space, two binary dihedral spaces and three spaces belonging to cyclic groups. In panel 
(a) the expectation values of the angular power spectra 5Tf are shown for the first 
three multipole moments 5T| (solid curve), 5T 3 2 (dashed curve) and 5T% (dotted curve). 
(Panel (b) in figures [T] to |H1 will be discussed below.) The corresponding values measured 
by the WMAP team are indicated as straight horizontal lines. With respect to the large 
scale anisotropy, the best values for f2 tot are obtained by choosing those values which 
yield the best agreement with the data. For example, consider the binary tetrahedral 
space shown in figure [T] Since the density Q tot should be as close to one as possible, 
one chooses the first minimum of <5T|, i. e. the range Q t ot = 1-06 . . . 1.07. In this range 
the value of 5T 3 2 is also near to the observed one. One should note that these values are 
only the expectation values such that for a given realization one has also to take into 
account the cosmic variance. However, as the simulations for such realizations show, 
the probability for a given spherical space increases significantly when the expectation 
values are already near to the observed values. In this way one gets for the binary 
octahedral space (figure EI) the range Q tot = 1.03... 1.04 and for the dodecahedron 
(figure H fitot = 1.015... 1.02. 

For none of the considered binary dihedral spaces is a comparably good agreement 
found. We have computed the CMB anisotropy for the binary dihedral groups D\ m 
with m = 2, 3, 4, 5, 10, 20, and 30. In figures 0] and El we show as two typical examples 
the cases m = 2 and m = 5, respectively. The first multipole moment 5T 2 2 decreases 
for all considered groups D\ m only for unrealistically high values of Q t ot > 1-1- In this 
parameter range the values for 5Tf are very large. Thus a binary dihedral topology 
seems to be very unprobable as a possibility for our Universe. 

The cyclic groups Z m are also much worse compared to the binary tetrahedral 
space, the binary octahedral space, or the dodecahedron. We have computed all groups 
Z m with m < 20 and a lot of examples up to m = 500. As three examples we present 
in figures El to |H1 the models for m — 1, 2, and 6, respectively. The group Z\ has only 
the identity as a group element and thus leads to the usual spherical space S 3 . One 
observes the known difficulty of the concordance model, i. e. too large values for the 
two lowest multipole moments. The next group shown is Z 2 leading to the projective 
space V 3 , also known as elliptic space, which has historically played a special role as an 
example of an alternative to the spherical space S 3 . For Q tot < 1.1 the first multipole 
moments behave very similar to the former case with m = 1 such that this topology is 
not a better match than the concordance model. Models with larger groups Z m lead to 
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a quadrupole moment which increases with increasing Q to t, in general. As an example, 
figure |H1 displays the result for Z 6 . The quadrupole moment <5T| shows a very small 
minimum around Q to t — 1-05 and increases for higher values of f2tot- Although there are 
minima in many cyclic models, they never suppress the power of 5T 2 2 even to the level 
of the observed value for 5T%, see figure El The fact that lens spaces do not fit well the 
CMB anisotropy is already discussed in [3U] where this behaviour is ascribed to the non 
well-proportioned fundamental cells. 

Therefore, from all spherical spaces, only the binary tetrahedral space, the binary 
octahedral space, and the dodecahedron display for the given parameter ranges the 
observed suppression in the large anisotropy power. 

Up to now, we have only discussed the angular power spectrum 5Tf. Let us 
now come to the temperature two-point correlation function C(i9), which is defined as 
C{d) := (5T(n)5T(n')) with n ■ n' = cos$. It can be computed from the multipole 
moments (jl4j) under the assumption of statistical isotropy as 

C(tf) ~ — £(2Z + l)C,^(cos0) . (25) 

1=2 

This quantity is well suited in order to measure the large scale power as emphasised 
in jS] where the observations are compared with the theoretical models. A comparison 
of C($) observed by WMAP with the concordance model can also be found in figure 
1 of reference (In we have derived an analytic expression for C sw ($), i.e. 

for the SW contribution in the case of the dodecahedron, which after multiplication by 
N/120 holds for all homogeneous spherical space forms.) The correlation function C{d) 
displays a surprisingly low CMB anisotropy on large angular scales d > p, which can be 
quantified by the S(p) statistic [TU] 

/cos p 
\C{d)\ 2 dcosd (26) 

which is discussed for the first-year WMAP data in [H] for p = 60°, and it is found 
that only 0.3% of the simulations based on the concordance model ri have lower values 
of S^O ) than the observed value 5(60°) = 1644. The S(p) statistic is shown for 
the above discussed spherical spaces in figures [T] to |S1 in panel (b) for p = 60° (solid 
curves) and p = 20° (dotted curves). The corresponding WMAP values are indicated as 
straight horizontal lines. The inspection of these figures leads to the same result as the 
above discussion of the first three angular power moments 5Tf. The binary tetrahedral 
space possesses a sufficiently low power in the range Q tot = 1.06... 1.07, the binary 
octahedral space in the range Q to t = 1-03. . . 1.04, and the dodecahedral space in the 
range Qtot — 1-015 . . . 1.02. For the other spherical spaces no comparable agreement is 
found as can be seen in the case of the two binary dihedral spaces (figures EJb) and 
Efb)) and the three cyclic groups (figures E^b) to|Bfb)). 

Now, we would like to demonstrate the strong influence of the radial function 
R%(t) on the suppression of power in the case of the quadrupole moment I = 2 already 
discussed above. The quadrupole suppression gets stronger, the higher the value of the 
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1.0 1.05 1.1 1.15 1.2 !l tot 1.0 1.05 1.1 1.15 1.2 Q, 

Figure 9. Panel (a) shows the dependence of the conformal distance tsls to the 
surface of last scattering on the density Sltot ■ Panel (b) shows the square of the radial 
function i?^ ; (r) for the first eigenvalue occurring in (|24|l for the quadrupole moment 
I = 2 for the binary tetrahedral space ((3 = 7), the binary octahedral space (,9 = 9), 
and the dodecahedron (/? = 13). 

first contributing eigenvalue is, see equation (fT7j) . In addition, it is seen from equation 
(fTTj) that the term belonging to the first eigenvalue ft is multiplied by R%(t). Thus the 
contribution of the first eigenvalue is eliminated for those values of f2 tot which belong to 
a value of r SLS at which -R/3z(t SL s) is zero. The dependence of r SLS on the density fl tot 
is shown in figure Efa) for our choice of cosmological parameters and is well described 
by t sls = 3.32-\/fitot _ 1 f° r 1 < ^tot < 1-1. In figure Efb) the square of the radial 
function -R^(t) is shown for the first eigenvalue of the binary tetrahedral space (ft = 7), 
the binary octahedral space (ft — 9), and the dodecahedron (ft = 13) for I = 2. One 
observes that the quadrupole suppression due to the first zero of the radial function is 
maximal in the case of the binary tetrahedral space at Q t ot — 1-064, at f2 tot ~ 1.038 
for the binary octahedral space, and at Q tot ~ 1.018 for the dodecahedral space. This 
matches perfectly well to the previously found intervals on which these models show 
a strong anisotropy suppression. Thus, the radial function has an important influence 
on the suppression for a given spherical topology. Note that the zero of the radial 
function eliminates many eigenfunctions due to the high degeneracy, e. g. the first 13 
eigenfunctions in the case of the dodecahedron. This does not happen in such a dramatic 
way in models with negative curvature, where one has no degeneracies at all, i. e. all 
eigenvalues have multiplicity one, in general. Then the radial function can only suppress 
a single eigenfunction and not a "cluster" of them. 

The angular power spectrum <5Tp as well as the S(p) statistic lead to the conclusion 
that there are three best candidates with respect to spherical spaces, i. e. the binary 
polyhedral spaces S 3 /T*, S 3 / O*, and S 3 / 1*. In figureslTUltoll2lwe show the temperature 
fluctuation ST/T in the Mollweide projection for these three topological spaces, where 
exactly those values of f2 to t are used which lead to a strong suppression of large scale 
power. In these calculations we have used for the three spaces all modes below the wave 
number cut-off /3 max = 155, 161 and 185, respectively. 
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Figure 10. The temperature fluctuation ST/T of one realization for the binary 
tetrahedral group T* is shown (/3 m ax = 155). The cosmological parameters fitot = 
1.065, fl A = 0.785 and h = 70 are used. 




Figure 11. The temperature fluctuation 5T/T of one realization for the binary 
octahedral group O* is shown (/3 m ax = 161). The cosmological parameters fitot = 
1.038, Q A = 0.758 and h = 70 are used. 

In figures El to ED we show the angular power spectrum 5Tf for the binary 
polyhedral spaces S 3 /T*, S 3 /0*, and S 3 /I*, where the same cosmological parameters 
as in figures HO to H3 are used. The la deviations are computed along the lines of 
|15j . Since the distributions for the lowest multipole moments are asymmetric, i.e. not 
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Figure 12. The temperature fluctuation ST/T of one realization for the binary 
icosahedral group I* is shown (/3 m ax = 185). The cosmological parameters fitot = 
1.018, n A = 0.738 and h = 70 are used. 
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Figure 13. The angular power spectrum 5Tf is shown for the binary tetrahedral 
group T* (open circles) using the same cosmological parameters as in figure ^| The 
angular power spectrum is shifted by Al — 0.25 in order to enable a comparison with 
the first-year WMAP data (full diamonds). The la errors are shown. 



Gaussian, these error bars should only be considered as providing the order of magnitude 
of fluctuations in individual realizations. In order to facilitate a comparison with the 
WMAP data, shown as full diamonds together with their la errors not including the 
cosmic variance, the spectra of the binary polyhedral spaces are shifted by Al = 0.25. 
The angular power spectra 8Tf for these three binary polyhedral spaces are very similar 
such that one is faced with a topological degeneracy with respect to 5Tf. All three 



CMB Anisotropy of Spherical Spaces 



18 



2000 




25 I 30 



Figure 14. The angular power spectrum 5Tf is shown for the binary octahedral 
group O* (open circles) using the same cosmological parameters as in figure ITTI 




25 I 30 



Figure 15. The angular power spectrum STj 2 is shown for the binary icosahedral 
group /* (open circles) using the same cosmological parameters as in figure FPU 



spectra display a good agreement with the WMAP data. 



4. Conclusion 



In this paper we analyse the CMB anisotropy of homogeneous 3-spaces of constant 
positive curvature which are multi-connected and are given by the quotient of S 3 by a 
group T of covering transformations, i.e. M. = S 3 /T. The motivation is provided by 
the surprisingly low power in the CMB anisotropy at the largest scales as measured by 
COBE and WMAP and the fact that the mean value of f2 to t reported by WMAP is 
1.020 which hints to a positively curved Universe. In order to explain this low power, 
one could modify the primordial power spectrum P$(/3), e.g. by carefully choosing the 
inflationary scalar potential, or by resorting to multi-connected space forms which give 
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a low CMB anisotropy at the largest scales due to missing modes compared to the 
simply-connected S 3 , in general. 

We study all types of homogeneous multi-connected spherical space forms and find 
no agreement for the cyclic groups Z m which show an enhanced power at the largest 
scales despite their small volumes. Also the binary dihedral groups D\ m do not lead to 
models with a suppression significantly stronger than the simply-connected S 3 universe. 
Thus these models do not seem to provide viable space forms as a model for our 
Universe. This contrasts to the remaining three space forms, the binary tetrahedral, 
the binary octahedral as well as the dodecahedral space forms which show a sufficiently 
strong suppression of large scale power compared to the simply-connected S 3 universe. 
The binary tetrahedral space requires a density O tot in the range 1.06... 1.07. Since 
the WMAP team reported Qtot = 1-02 ± 0.02 this model is probably in conflict with 
the observations. For the two remaining models the density Qtot should be in the 
range Q tot = 1.03... 1.04 for the binary octahedral space, and f2 tot = 1.015... 1.02 
for the dodecahedral space. These values are compatible with the current observations. 
Furthermore, we would like to remark that the binary octahedral space displays a slightly 
stronger suppression of power than the dodecahedral space, as a comparison of figures 
121 and figures El reveals. 

A unique signal for a particular topology is provided by the so-called circles-in-the- 
sky-signature proposed in [31]. Along two circles on the sky which are mapped onto each 
other by the group T, the ordinary Sachs- Wolfe effect produces the same temperature 
signal. If there would be no Doppler and integrated Sachs-Wolfe contribution, see 
equation (J2H), which disturb this signal, one would expect a clear sign for a given 
topology if present. In 14 we study the influence of the latter two contributions on the 
circles-in-the-sky-signature and find that the degradation of the signal is strong enough 
such that the topology signal can be swamped. Therefore, the fact that in [32] no circles 
are found in the WMAP sky maps, does not necessarily exclude the binary octahedral 
or the dodecahedral space as viable models for our Universe. In a forthcoming paper 
we will study the circles-in-the-sky for the three best multi-connected spherical space 
forms, and in particular shall discuss whether a combined circle search on all circles 
simultaneously of a given topology can overcome the degradations. 
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